Tensor interaction contributions to single-particle energies 
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We calculate the contribution of the nucleon-nucleon tensor interaction to single-particle energies 
with finite-range G matrix potentials and with zero-range Skyrme potentials. The Skx Skyrme 
parameters including the zero-range tensor terms with strengths calibrated to the finite-range results 
are refitted to nuclear properties. The fit allows the zero-range proton-neutron tensor interaction 
as calibrated to the finite-range potential results and that gives the observed change in the single- 
particle gap e(hn/ 2 )-e(g7/2) going from 114 Sn to 132 Sn. However, the experimental I dependence of 
the spin-orbit splittings in 132 Sn and 208 Pb is not well described when the tensor is added, due to 
a change in the radial dependence of the total spin-orbit potential. The gap shift and a good fit to 
the ^-dependence can be recovered when the like-particle tensor interaction is opposite in sign to 
that required for the G matrix. 
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The tensor force between nucleons is important for the 
single-particle energies spacing and the shell structure of 
nuclei obtained from shell-model configuration interac- 
tion models of nuclei [l|, @- But, except for an early 
exploratory work Q , its role in Hartree-Fock models has 
been neglected. Results for relative energy shifts have 
recently been obtained with a finite-range tensor in the 
Gogny model [J]. In this letter we present the first sys- 
tematic results for a Skyrme-type interaction with a zero- 
range tensor interaction from a global fit to nuclear data 
including data for single-particle energies. We start by 
calibrating the strength of the zero-range tensor with re- 
sults obtained from a G matrix interaction. Then the 
Skyrme parameters plus the tensor parameters are var- 
ied to obtain a best fit to data. We find that the fit 
will allow an isovector tensor strength similar to that ex- 
pected from the G matrix. But the data set prefers an 
isoscalar tensor that is much smaller than expected from 
the G matrix. The reason is traced to a difference in be- 
tween the Skyrme spin-orbit and tensor radial forms in 
the doubly closed-shell nuclei 132 Sn and 208 Pb. 

First we calculated the contribution to the single- 
particle proton energies for single-particle states above 
the Z = 50 closed shell from the tensor part of the 
Hosaka-Kubo-Toki (HKT) G matrix @. HKT is a one- 
boson exchange potential that reproduces the G matrix 
elements obtained from the Paris potential. This tensor 
interaction has the form: 
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Xi = r/ri where are the range parameters. This con- 
sists of the one-pion exchange potential with r„ = 1.414 
fm, [W«,T=o/W n , T =i] = -3, and W„, T =i = 3.49 MeV, 
plus short-range potential with r s = 0.25 fm and with 
the strengths W determined from the G matrix elements: 
W s , T =o = 3105 MeV and W s , T =i = -1382 MeV. 

The contributions to the proton single-particle states 
for 132 Sn are shown in Table I. They were obtained with 
harmonic oscillator radial wavefunctions with Huj=7.87 
MeV. (The results with the tensor part of the M3Y po- 
tential [(| are the same as HKT within about 3%). The 
contribution to the single-particle energy of the valence 
proton in orbital k — (n,l,j) from the core protons is 
obtained from: 



(2j + l)Ei p = J2(2J+lKk> 

k',J 

and from the core neutrons from: 



J,x=u 



(2) 



(2j + l)E{ n = E 1± | M1 (2^ + 1)W (3) 

k',J,T 

where the two-body matrix elements are 
Vlk',j,T=<k,k',J,T\V t | k,k',J,T>. 

In the sum over the core orbitals k' the contributions 
from the sum of j'C, = f + 1/2 and j< = £' - 1/2 or- 
bital pairs cancel when both are filled as shown in Eq. 
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TABLE I: Contributions of the tensor finite-range G matrix 
interaction to single-particle proton energies in 132 Sn. The re- 
sults are given for E kp : contribution from the Ogg/2 proton or- 
bital, i?£ n (100): contribution from the Ogg/2 neutron orbital, 
Sfe n (114): i?J, n (100) plus the contribution from the O517/2 and 
14/2 neutron orbitals, E\ n : _EJj. n (114) plus the contribution 
from the 14/2 and O/in/2 neutron orbitals. The short-range 
tensor contribution is given by the "s-only" results. 
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k = (ntj) 


Ek P 


sL(ioo) 


EL(1U) 


Ekn 


Ek p 






(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


7T + S 


0<?7/2 


-0.458 


-1.009 


-0.135 


-1.032 


-1.490 




14/2 


0.078 


0.180 


0.395 


0.218 


0.296 




14/2 


-0.118 


-0.270 


-0.593 


-0.328 


-0.446 




Ohu/2 


0.308 


0.688 


0.109 


0.848 


1.156 


s only 


0ff7/2 


0.251 


0.408 


0.072 


0.465 


0.716 




14/2 


-0.060 


-0.097 


-0.162 


-0.100 


-0.160 




14/2 


0.090 


0.145 


0.243 


0.150 


0.240 




O/111/2 


-0.191 


-0.310 


-0.050 


-0.397 


-0.588 


7T + S 
S 


0ff7/2 


-1.82 


-2.48 


-1.88 


-2.22 


-2.08 




14/2 


-1.31 


-1.86 


-2.44 


-2.18 


-1.85 




14/2 


-1.31 


-1.86 


-2.44 


-2.18 


-1.88 




0/lll/2 


-1.62 


-2.22 


-2.18 


-2.14 


-1.97 



(4) of [2]. Thus the E l kq are zero for LS closed cores. For 
non-LiS closed cores with a pair of valence orbits with 
j > = 1 + 1/2 and j< =1—1/2 the energy shifts are 
(2j> + l)El q = -(2j < + l)E~l <q , which means that for 
a given £ value the tensor interaction with the core con- 
tributes to the effective spin-orbit splitting (see Eq. (4) 
of 0). The short-range contribution (s) to the energy 
shifts are given in the middle part of Table I. One ob- 
serves the change in sign that is related to the partial 
cancellation of the 7r-exchange potential by the short- 
range potential. The variation of the ratio over several 
valence orbitals shown at the bottom of Table I is a mea- 
sure of how well the finite-range tensor can be approxi- 
mated by a zero-range form. This is important since the 
zero-range approximation leads to an analytic form for 
the tensor density functional, and an efficient implemen- 
tation in the Skyrme Hartree-Fock method [|[ . The ratio 
varies depending on which orbits are filled in the core, 
up to a factor of two. But in the case of total energy 
for protons in 132 Sn the orbit dependence in the ratio is 
small. 

One observes from Table I that the tensor interaction 
results in a change in the O57/2-O/111/2 gap going from 
114 Sn (e.g. where only the 14/2 an d 0<77/ 2 orbitals are 
filled) to 132 Sn of 1.64 MeV. 132 Sn is one of the best 
doubly magic nuclei, and the lowest levels in 133 Sb are 
thus taken as single-particle states for adding a proton to 
132 Sn, although the experimental measurement of spec- 
troscopic strength from one-proton transfer reactions has 
not yet been carried out. The single-particle proton en- 



TABLE II: Proton single-particle energies in Sn and Sn. 
The gap is the energy difference between Og/7/2 and O/in/2. 



Nucleus 




exp 


exp 


OKX 


oKxta 


OKXtD 






lowest J 


centroid 












(MeV ) 


(MeV ) 


(MeV J 


(MeV J 


(MeV ) 


" 2 Sn 


0^7/2 


-9.68 


-9.68 


-9.87 


-10.82 


-9.86 




14/2 


-8.72 


-8.72 


-9.20 


-9.22 


-9.30 




14/2 


-6.97 


-6.97 


-7.38 


-7.50 


-7.14 




2S1/2 






-6.82 


-6.93 


-6.78 




04l/2 


-6.89 


-6.89 


-6.92 


-6.08 


-6.66 




gap 


2.79 


2.79 


2.84 


4.52 


3.20 


114 Sn 


Off 7/ 2 


-3.01 


-2.48 


-2.78 


-2.70 


-1.59 




14/2 


-3.73 


-3.02 


-2.91 


-2.81 


-2.66 




14/2 


-2.66 


-1.63 


-0.88 


-1.06 


-1.16 




2S1/2 


-2.96 


-1.32 


-0.81 


-0.87 


-0.85 




0/lll/2 


-2.43 


-1.45 


0.19 


0.11 


-0.76 




gap 


0.58 


1.03 


2.97 


2.81 


0.83 



ergies for 132 Sn are given in Table II. 

Proton transfer experiments have been carried out for 
114 Sn to 115 Sb. But the interpretation of the experimen- 
tal results in terms of single-particle energies is not so 
simple since the neutron configuration in 114 Sn is not 
magic with significant configuration mixing between the 
lowest neutron orbits of 0(77/2 and 14/2 and the upper 
orbits of 14/2> 2s 1 / 2 and 0/in/ 2 . To estimate the ef- 
fect of splitting of single particle strength we carry out 
a large-basis shell-model calculations that includes up to 
three neutrons being excited from the lower to the upper 
orbits with the renormalized G matrix interaction from 
[7j. The spectroscopic strength obtained is shown in fig- 
ure 1. One observes the lowest J = j states contains the 
largest fraction of single-particle strength, but that there 
is significant spreading to higher energy. The isolation 
of one large part of the spectroscopic strength into the 
lowest state is consistent with experimental observation 
Q. This spreading is due to coupling with the neutron 
vibrations within the 0177/2, Id, 2s, hn/2 model space as 
well as isospin splitting of the strength to the T> = 15/2 
states. The centroid energies obtained from the results of 
Fig. 1 are within about 10 keV of the single-particle en- 
ergies obtained with the simplest [i/,0gy 2 ,ld®/ 2 )[ir nt ij] 
configuration. This simple configuration is the one as- 
sumed for the finite-range tensor contribution discussed 
above as well as for the Skyrme Hartree-Fock calculations 
to be discussed below. (This configuration does not have 
good isospin, but configuration mixing restores isospin.) 
In order to estimate the proton single-particle energies in 
115 Sb we add on to the separation energy of the lowest 
states of a given J = j in the experimental spectrum a 
correction based on the configuration mixing results of 
Fig. 1, giving the experimental centroid energies in col- 
umn 3 of Table II. 
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FIG. 1: Single-particle proton spectroscopic factors for 114 Sn 
to 115 Sb from the shell-model calculations. The lines for each 
n, £, j value correspond to the cumulative sum of spectroscopic 
strength as a function of the energy-difference Ef — Ei. The 
centroid energies are indicated by the dashed lines. 

In addition to configuration mixing within the 
0<?7/2> Id, 2s, hn/2 model space, one could consider the 
effects of protons excited across the Z — 50 shell gap 
and neutrons excited across the N — 50 and N = 82 
shell gaps. The proton excitations have a direct affect 
on enhancing the B(E2) values for low-lying 2 + states of 
Sn @. This type of configuration mixing affects single- 
particle energies for all of the Sn isotopes including those 
for 133 Sb. However, this is already partly accounted for 
in Skx with its empirical effective mass of near unity [Io| . 
The enhancement of the bare G matrix effective mass of 
m*/m = 0.6 — 0.7 towards its empirical value of near 
unity in Skx can be attributed to coupling with the mul- 
tipole vibrations of the core [ll|, 

The experimental centroid energies are compared with 
the results of Skyrme Hartree-Fock calculations. The Skx 
interaction [To| was obtained from a fit to binding en- 
ergies, rms charge radii and single-particle energies in- 
cluding those for the 132 Sn core [lCj]- Thus the rather 
good agreement between the experimental and Skx pro- 
ton single-particle energies for 132 Sn is not an accident. 
The 114 Sn data were not included in the original Skx fit 
since, as discussed, this does not have a magic neutron 



number. However, the agreement with the centroid en- 
ergies is not bad except for the 0hn/ 2 . 

The gap between the O/in/2 and 0g 7 / 2 single particle 
energies is particularly sensitive to the tensor interaction. 
Experimentally the gap changes by 1.76 MeV; from 2.79 
MeV in 132 Sn to 1.03 MeV in 114 Sn. The gap for Skx is 
about the same for 114 Sn and 132 Sn and this is similar 
to what we find for other Skyrme interactions. But the 
finite-range tensor interaction discussed above leads to 
about a gap shift of 1.64 MeV - close to the observed 
shift of 1.76 MeV and to the values shown in Fig. 4(d) 
of and Fig. 4 of Thus we are motivated to add a 
tensor interaction to the Skyrme functional. 

We use the zero-range form of the tensor potential 
given by Stancu et al. Q (Eq. 1 of their paper). The 
zero-range tensor gives an additional contribution to the 
Skyrme spin-orbit term of the form: 

AW n = aJ n + (3J P , AW P = aJ p + /3J n (4) 

where the coefficients a and (3 (given in units of MeV fm 5 ) 
come from the zero-range form of the tensor interaction 
(at and f3t) as well as from the exchange part of the 
central interaction: 

otc = -t 2 ) - ^{hxi + t 2 x 2 ) (5) 

and 

C = --(t%xx+t 2 x 2 ). (6) 
o 

The J q are the spin densities defined by 

J i(r) = S(2ja+l)[7a(ja+l)-^(^+l)-|i^(r) 

a 

where the sum is over the occupied orbits with proton 
(q = p) or neutrons (q = n). 

We start with the Skx interaction [l(| and the data 
base that was used to determine its parameters. Skx 
is the only interaction for which a large number of ex- 
perimental single-particle energies where used to con- 
strain the parameters. As a baseline for our new fits, 
Skx gives a % 2 value of 0.60 when the parameters 
tp,ti,t 2 , t 3 , xq, xi,x 2l X3 and W are fitted to the data set 
of [l0| . For this original Skx the a c and (3 C terms were not 
included. If they are included, the % 2 increases slightly 
to 0.62 and the central-exchange values are a c = 24 and 
Pc = -23. 

As For the tensor contribution, the initial set of at and 
flt parameters were chosen to reproduce the calculated 
E pk and E^ k values for the 0*77/2 proton orbit from the 
finite-range G matrix given in Table I. The results are 
at=60 and /?t=110. These are larger than Skx central- 
exchange values of a c — 24 and j3 c — -23, but both should 
be considered for the total and the refits we carry out will 
include the effects of a c and f3 c . For comparison with 
Stancu et al. [3[ , our values of at and f3 t are close to the 



TABLE III: Contributions to proton single-particle energies 
in 132 Sn for the Skxta and Skx Hartree-Fock calculations. E* 
is the total tensor contribution with a=93 and /3=94, with 
the zero-range tensor interaction contribution for a t —60 and 
/3 t = 110 shown in brackets. E s ° is the spin-orbit potential 
contribution to the energy. The total for Skxta E\ v + E\ n + 
E so is compared with the spin-orbit for Skx. 











total 






Skxta 


Skxta 


Skxta 


Skxta 


Skx 




(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


O37/2 


-0.723 (-0.476) 


-0.828 (-0.984) 


3.20 


1.65 


2.64 


lds/2 


0.117 (-0.079) 


0.089 (0.094) 


-0.83 


-0.62 


-0.68 


ld 3 /2 


-0.181 (-0.121) 


-0.158 (-0.171) 


1.34 


1.00 


1.11 


Ohu/2 


0.640 (-0.421) 


0.864 (1.034) 


-3.84 


-2.34 


-3.19 



those the y es timate from the interaction of Sprung and 
Banerjee [l| with q = 1.0 far 1 (Table 1 of gjj). 

Then the Skyrme parameters to , t\ , ti , £3 , xq , x\ , x-i , £3 
and W were refit for these fixed values of at and fit with a 
resulting set of parameters called Skxta. The x 2 value in- 
creased significantly to 1.50. The contributions from the 
zero-range tensor, central-exchange and spin-orbit inter- 
actions to the proton single-particle energies in 132 Sn are 
shown in Table III. The central exchange values from the 
fit are ev c =33 and /3 C =-16. 

The single-particle energies for orbitals around 132 Sn 
obtained with Skx and Skxta are shown in Figs. 2 and 
3 respectively. Comparison of these figures shows that 
the x 2 increase is due to a poorer ^-dependence of the 
spin-orbit splitting for Skxta compared to Skx. This can 
be traced to a difference in the radial functional form 
of the spin-orbit contributions that are shown in Fig. 4. 
The tensor contribution peaks at a 0.5 fa smaller radius 
compared to the Skyrme spin-orbit potential. Since the 
tensor contribution with a t —60 and fj t —H0 is opposite in 
sign to the normal spin-orbit contribution, the strength 
of the Skyrme spin-orbit parameter W has to increase by 
about 20% to recover an overall fit to the single-particle 
energy data. But the ^-dependence of the experimental 
single-particle energies are better reproduced with the 
Skyrme spin-orbit shape. The quality of the Skx and 
Skxta fit results for single-particle energies around 208 Pb 
are similar to those we show for 132 Sn. 

Given that /3t=110 is needed to reproduce the shift of 
the 0g 7 / 2 -0/iii/2 gap going from 114 Sn to 132 Sn, we next 
fix /3 t = 1 1 and include at in the Skyrme fit to data. We 
recover a good fit of % 2 =0.63 with a t = —118 for a pa- 
rameter set we call Skxtb. In general we find a good 
fit with values of j3t in the range of to 110 as long as 
a t « — (3f This happens because the proton and neu- 
tron contributions then cancel in the jj closed shell nu- 
clei 132 Sn and 208 Pb giving the good reproduction of the 
single-particle energies from the Skyrme spin-orbit shape. 
The Skxtb single-particle energies are compared with ex- 




FIG. 2: Comparison of experimental and theoretical single- 
particle energies in 132 Sn for the Skx interaction. 




FIG. 3: Comparison of experimental and theoretical single- 
particle energies in 132 Sn for the Skxta interaction. 
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FIG. 4: Spin-orbit potentials for protons in 132 Sn for the 
Skxta interaction: the Skyrme spin-orbit potential (crosses); 
the zero-range tensor contributions from the core protons 
(dashed line) and core neutrons (dotted line); and the to- 
tal (solid line). The total also includes the smaller a c and /3 C 
terms from the exchange part of the central interaction. 

periment in Table II. Skxtb gives a best account of the 
132 Sn single-particle energies and the 0^7/2 — Ohn/2 gap 
shift. However, the absolute single-particles energies in 
114 Sn still differ from experiment by as much as one MeV. 
In conclusion, we find that the finite-range tensor inter- 



action is important for the Ogr/2 — Ohu/2 gap shift. How- 
ever, a zero-range implementation of the tensor interac- 
tion in the Skryme interaction is problematic. The radial 
form of the tensor contribution to the spin-orbit potential 
does not give a good reproduction of the ^-dependence of 
the spin-orbit splittings in 132 Sn and 208 Pb. Reproduc- 
tion of the observed 0(77/2 — Ohn/2 gap shift plus a good 
fit to absolute single-particle energies in 132 Sn and 208 Pb 
requires fit ~ 110 for the proton-neutron tensor interac- 
tion (consistent with the G matrix value), and at ~ —fit- 
for the T = 1 tensor interaction between like particles, 
which is opposite in sign to the G matrix value of at ~ 60. 
Although our finite-range calculations indicate that the 
zero-range approximation may be adequate, further in- 
vestigation is required. Also we need to understand the 
role of correlations (coupling to vibrations) and three- 
body forces on the effective tensor interactions in nu- 
clei. The central and part of the Skyrme functional also 
need to be constrained and extended to reproduce realis- 
tic properties of nuclear matter [3] ■ These changes may 
lead to different values for a c and fi c than those obtained 
with Skx that also need to be taken when the tensor in- 
teraction is included. 
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